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Abstract 

We review the present status of simulations of top production and decay in 
the HERWIG Monte Carlo event generator. We show the phenomenological 
impact of the recently-implemented matrix-element corrections to top decays 
for e~^e~ collisions at 360 GeV and discuss possible further improvements for 
studies of future experiments at the Linear Collider. 
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1 Introduction 



The study of the top quark phenomenology will be among the main topics of interest 
for the future experiments at the Linear Collider. In order to perform such analyses, it 
will be essential to have reliable Monte Carlo simulations of top production and decay. 
Standard event generators |jl|,|^ simulate parton showers in the soft /coUinear approxi- 
mation, leaving empty regions in the phase space ('dead zones') where, according to the 
exact matrix element, the radiation should be suppressed but not completely absent. 
Parton showers need therefore to be supplemented by matrix-element corrections. 

As far as the HERWIG parton shower is concerned, matrix-element corrections to 
e+e~ annihilation 0], deep inelastic scattering [|], top decays ^ and Drell-Yan processes 
have been implemented following the general prescriptions given in |]^. In section 
2 we shall review the HERWIG parton shower algorithm for top production and decay 
in e~^e~ collisions and discuss the recent implementation of matrix-element corrections 
to top decays. In section 3 we shall show some phenomenological results at the Linear 
Collider, for a centre-of-mass energy of 360 GeV. In section 4 we shall make some 
concluding remarks on the results and possible improvements of this work. 



2 The HERWIG parton shower algorithm for top 
production and decay in e+e~ coUisions 

For the emission of one more parton (gluon) in top production in e~^e~ annihilation 
{e~^e^ tig) and top decay {t bWg), the elementary probability implemented in the 
HERWIG event generator is the general result for parton radiation in the soft / collinear 
limit g: 

= (1) 

r 27r As{q^,qi) 

The shower is ordered according to the variable = E'^C,, E being the energy of the 
parton that splits and ^ = pi ■ P2/ {E1E2), where pi and p2 are the momenta of the 
two outgoing partons; z is the energy fraction of the emitted parton (gluon) relative to 
the incoming one. In the massless approximation, ^ = 1 — cos^, 6 being the emission 
angle, in such a way that ordering according to g^, in the soft limit, results in angular 
ordering. As{q'^,q'^) is the Sudakov form factor, expressing the probability that no 
resolvable radiation is emitted from a parton whose upper limit on emission is g^, with 

being a cutoff on transverse momentum. The ratio of Sudakov form factors in (|l|) 
sums up all virtual contributions and unresolved emissions. 

The parton shower variables z and ^ are frame-dependent, nevertheless it is possi- 
ble to prove that colour coherence implies that the values gmax of any pair of colour- 
connected partons are related via gimaxgjmax = Pi -Pj, which is Lorentz- invariant. For 
top production, as for most of the HERWIG processes, symmetric limits are chosen, 
i-e. gimax = 9jmax = y^Pi"^ and the energy of the parton which initiates the cascade is 
equal to Ei = g^max- Ordering according to g^ leads to the condition ^ < 1. The phase 
space region 1 < ^ < 2 is not populated at all by the parton shower algorithm. In M 
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matrix-element corrections are implemented to e+e^ —>■ qq processes: the 'dead zone' of 
the phase space is populated according to a probability distribution obtained from the 
exact 0{as) matrix-element calculation of the process e~^e~ qqg ('hard correction') 
and the emission in the already-populated region is corrected using the exact amplitude 
for every emission that is the hardest so far ('soft correction'). 

Top decays are treated in a somewhat different way in HERWIG, since the top quark 
rest frame is chosen to perform such a decay This means that Et = nit = Itmax and 
Eb = Pt ■ Pb/^t- Being at rest, the top quark is not allowed to radiate gluons in the 
decay stage, while the b quark emits soft gluons for ^ < 1, i.e. < 9 < n/2. The result 
is that the soft phase space is not completely populated by HERWIG as all the soft 
gluons which should be emitted in the backward hemisphere 7i/2 < 6 < n are missed. 
For e~^e~ annihilation we have two back-to-back partons, q and q, that are capable of 
emitting soft gluons for ^ < 1, in such a way that each of them radiates in the region 
the other is missing out and the whole of the soft phase space is filled. In it is shown 
that, though neglecting the soft radiation in the backward hemisphere, the total energy 
loss due to gluon radiation is roughly right. However, we may have serious problems 
when dealing with angular differential distributions. 

For the reasons previously mentioned, the implementation of matrix-element correc- 
tions to top decays is not a straightforward extension of the method applied for the 
purpose of other processes. In order to avoid the soft singularity, we set a cutoff on the 
energy of backward gluons, whose default value is -Emin = 2 GeV. In [^j we calculated the 
phase space limits for the decay t{q) — ^ b{P2)W{pi)g{p3) with respect to the variables 

2p2 ■ P3 , . 

Xl = 1 (2) 

X3 = (3) 

In Fig. ID we plot the total and the HERWIG phase space. The differential width reads: 



2 



log Ya 

Figure 2: distributions according to HERWIG 6.1 (dashed line), HERWIG 6.0 (dot- 
ted) and according to the exact 0{as) calculation (sohd). 

I (PV 1 as ^ I { , 

r„d^, = (1 - a) (1 + i - 2a) ' ^ C'.^j|^|,..,-..,l,.-.,..-(.-„,(.-..,-.i] 

+ {1+^)X3(X1+X3-1)2+24(1-Xi)| , (4) 

where Tq is the width of the Born process t — >■ bW and a = m^/m^. The integral of the 
differential width over the 'dead zone' is divergent, because of the soft divergence. We 
use the differential width (^) to generate events in the missing backward phase space and 
in the already-populated forward hemisphere for every hardest-so-far emission capable 
of being the hardest so far. Both hard and soft corrections should be applied only for 
gluon energies larger than the cutoff value. 



3 Phenomenological results at the Linear Collider 

In order to test the reliability of the HERWIG algorithm provided with matrix-element 
corrections we wish to compare its phenomenological results with the ones obtained by 
calculating the exact first-order matrix element of the process: 

e+e- ^tt^ {bW+){bW-)g. (5) 

It is interesting to perform such an analysis at the centre-of-mass energy ^/s = 360 GeV, 
slightly above the threshold for tt production, in such a way that all the gluon emission is 
associated to top decays, the available phase space for radiation in the production stage 
being too small. In the last public version HERWIG 5.9 some bugs were found in the 
treatment of top decays, therefore, if we wish to investigate the effect of the implemented 
matrix-element corrections, it is better to compare the new version HERWIG 6.1 to an 
intermediate version, which we call HERWIG 6.0, where all the bugs are supposed to 
be fixed. 
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Figure 3: Invariant opening angle distributions according to HERWIG 6.1 (dashed line), 
HERWIG 6.0 (dotted line) and according to the exact 0(0:5) matrix element (solid). 



We study three-jet events according to the Durham algorithm [|TU| and plot the 
differential distributions with respect to the threshold variable and the minimum 
invariant opening angle AR = ^ {Acpy + (Ar])^ using HERWIG 6.1 and 6.0 and the 
exact 0{as) for a value as = 0.145 of the strong coupling constant. We also set the 
cuts Et > 10 GeV and AR > 0.7 on the resolved jets. Such an analysis was already 



performed in [|1^] and serious discrepancies were found when comparing the exact results 
with the ones obtained by the using of HERWIG before the implementation of matrix- 
element corrections to top decays. 

In Figs. and ^ we observe a remarkable impact of matrix-element corrections: 
HERWIG 6.1 generates more events at large with respect to the 6.0 version, which 
can be explained as due to the hard corrections which allow hard and large-angle gluon 
radiation. We also get a suppression at small ys, which compensates the enhancement at 
large 7/3, the total amount of radiation being roughly right even before matrix-element 
corrections. The soft corrections, applying the matrix-element distribution instead of 
the parton shower one in the already-populated forward region for every hardest-so- 
far emission, is responsible of the suppression at small y^. The agreement between 
HERWIG 6.1 and the exact calculation is pretty good at large 7/3, as it should be 
since the fixed-order calculation is reliable only if we are far from the soft and collinear 
divergences, which correspond to small y^, where, on the contrary, parton showers are 
more trustworthy. Similar comments hold for the AR plots as well. In |^ we also showed 
that the dependence on the chosen value of the cutoff i^^min is negligible after we apply 
the experimental cuts. 



4 Conclusions 

We have reviewed the HERWIG parton shower model for top production and decay in 
e^e~ annihilation and the main features of the method of matrix-element corrections, 
recently implemented to the simulation of top decays. We studied three-jet events at 
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a/s = 360 GeV using HERWIG before and after matrix-element corrections to top 
decays, and the exact first-order calculation. The results showed a marked effect of 
the improvement to the treatment of top decays and a good agreement with the 0{as) 
calculation in the phase space region where such an agreement is to be expected. We 
therefore feel confident that the new version HERWIG 6.1 should be a reliable Monte 
Carlo event generator to simulate top quark decay at the future Linear Collider. 

For the sake of completeness, we have however to say that the implementation of 
matrix-element corrections to the process e~^e~ qqg still needs some improvement 
since, in the determination of the phase space limits and in the application of the 
soft corrections, the qq pair is treated as if it was massless, which might not be a good 
approximation for the purpose of the top quark. We have been working on fully including 
these mass effects. Furthermore, it will be very interesting to investigate the impact of 
matrix-element corrections to top decays on the top mass reconstruction at the Linear 
Collider. This is in progress as well. 
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